Device Model for Graphene Bilayer Field-Effect Transistor 
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We present an analytical device model for a graphene bilayer field-effect transistor (GBL-FET) 
with a graphene bilayer as a channel, and with back and top gates. The model accounts for the 
dependences of the electron and hole Fermi energies as well as energy gap in different sections of 
the channel on the bias back-gate and top-gate voltages. Using this model, we calculate the dc and 
ac source-drain currents and the transconductance of GBL-FETs with both ballistic and collision 
dominated electron transport as functions of structural parameters, the bias back-gate and top-gate 
voltages, and the signal frequency. It is shown that there are two threshold voltages, Vth.i and 
V t h,2, so that the dc current versus the top-gate voltage relation markedly changes depending on 
whether the section of the channel beneath the top gate (gated section) is filled with electrons, de- 
pleted, or filled with holes. The electron scattering leads to a decrease in the dc and ac currents and 
transconductances, whereas it weakly affects the threshold frequency. As demonstrated, the tran- 
sient recharging of the gated section by holes can pronouncedly influence the ac transconductance 
resulting in its nonmonotonic frequency dependence with a maximum at fairly high frequencies. 

PACS numbers: 73.50.Pz, 73.63.-b, 81.05.Uw 



I. INTRODUCTION 

The features of the electron and hole energy spectra in 
graphene provide the exceptional properties of graphene- 
based heterostructures and devices 0, 0, [H, 13, HI, H| • 

However, due to the gapless energy spectrum, the 
interband tunneling [7| can substantially deteriorate 
the performance of graphene field-effect transistors (G- 
FETs) with realistic device structures @, ©, 01, ED]. To 
avoid drawbacks of the characteristics of G-FETs based 
on graphene monolayer with zero energy gap, the pat- 
terned graphene (with an array of graphene nanorib- 
bons) and the graphene bilayers can be used in graphene 
nanoribbon FETs (GNR-FETs) and in graphene bilayer 
FETs (GBL-FETs), respectively. The source-drain cur- 
rent in GNR-FETs and GBL-FETs, as in the standard 
FETs, depends on the gate voltages. The positively bi- 
ased back gate provides the formation of the electron 
channels, whereas the negative bias voltage between the 
top gate and the channels results in forming a poten- 
tial barrier for electrons which controls the current. By 
properly choosing the width of the nanoribbons, one can 
fabricate graphene structures with a relatively wide band 
gap [3 (see also Refs. [3E1H0J]) Recently, the de- 
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vice dc and ac characteristics of GNR-FETs were assessed 
using both numerical [lj] and analytical [17lllcill9| mod- 
els. The effect of the transverse electric field (to the 
GBL plane) on the energy spectrum of GBLs [13, EL H2] 
can also be used to manipulate and optimize the GBL- 
FET characteristics. A significant feature of GBL-FETs 
is that under the effect of the transverse electric field not 
only the density of the two-dimensional electron gas in 
the GBL varies, but the energy gap between the GBL 
valence and conduction bands appears. This effect can 
markedly influence the GBL-FET characteristics. The 
structure of a GBL-FET is shown in Fig. 1. In this pa- 
per, we present a simple analytical device model for a 
GBL-FET, obtain the device dc and ac characteristics, 
and compare these characteristics with those of GNR- 
FETs. 

The paper organized as follows. In Sec. II, we consider 
the GBL-FET band diagrams at different bias voltages 
and estimate the energy gaps and the Fermi energy in dif- 
ferent sections of the device. Section III deals with the 
Boltzmann kinetic equation which governs the electron 
transport at dc and ac voltages and the solutions of this 
equation. The cases of the ballistic and collision domi- 
nated electron transport are considered. In Sec. IV and 
Sec. V, the dc transconductance and the ac frequency- 
dependent transconductance are calculated using the re- 
sults of Sec. III. Section VI deals with the demonstration 
and analysis of the main obtained results, numerical esti- 
mates, and comparison of the GBL-FET properties with 
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FIG. 1: Schematic view of the GBL-FET structure 



those of GNR-FETs. In Sec. VII, we draw the main con- 
clusions. In Appendix, some intermediate calculations 
related to the dynamic recharging of the gated section 
by holes due to the interband tunneling are singled out. 



II. GBL-FET ENERGY BAND DIAGRAMS 

We assume that the bias back-gate voltage V& > 0, 
while the bias top-gate voltage Vt < 0. The electric po- 
tential of the channel at the source and drain contacts 
are <p = and tp = Vd, respectively where Vd is the 
bias drain voltage. The former results in the formation 
of a 2DEG in the GBL. The distribution of the electron 
density £ along the GBL is generally nonuniform due to 
the negatively biased top gate forming the barrier region 
beneath this gate. Simultaneously, the energy gap E g is 
also a function of the coordinate x (its axis is directed 
in the GBL plane from the source contact to the drain 
contact) being different in the source, top-gate, and drain 
sections of the channel (see Fig. 2). Since the net top- 
gate voltage apart from the bias component Vt comprises 
the ac signal component 6V(t), the height of the barrier 
for electrons entering the section of the channel under 
the top gate (gated section) from the source side can be 
presented as 



A(t) = A + <JA(t). 



(1) 



Depending on the Fermi energy in the extreme sec- 
tions of the channel, in particular, on its value, ef, in 
the source section and on the height of the barrier in this 
section Ao, there are three situations. The pertinent the 
GBL-FET energy band diagrams are demonstrated in 
Fig. 2. The spatial distributions of electrons and holes in 
the GBL channel are different depending on the relation- 
ship between the top-gate voltage Vt and two threshold 
voltages, V t h,i and V t h,2- These threshold voltages are 
determined in the following. 

When Vth,2 < Vth,x < Vt, the top of the conduc- 
tion band in the gated section is below the Fermi level 
(Fig. 2a). In this case, an n + -n-n + structure is formed in 
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FIG. 2: Band diagrams at different top gate bias voltages 
(V b > 0, V d = 0): (a) V th , 2 < Vth,i < V, (b) Vifc.a < V < 
Vth,i (depleted gated section), and (c) V t < Vtk,2 < Vth.i 
(gated section filled with holes), Panel (d) corresponds to 
V th ,2 <V t < V th ,i but with V d > 0. 



the GBL channel. At V t h,2 < V t < Vth.i, the Fermi level 
is between the top of the conduction band and the bot- 
tom of the valence band in this section (Fig. 2b). This 
top-gate voltage range corresponds to the formation of 
an n + -i-n + structure. If Vt < Vth,2 < Vth.i, both band 
edges are above the Fermi level (Fig. 2c), so that n + - 
p and p-n + junctions are formed beneath the edges of 
the top gate. In the first and third ranges of the top 
gate voltage ( "a" and "b" ranges) , the electron and hole 
populations of the gated section are essential. In the sec- 
ond range (range "b"), the gated section is depleted. In 
the voltage range "a" , the source-drain current is associ- 
ated with a hydrodynamical electron flow (due to effec- 
tive electron-electron scattering) in the gated section. In 
this case, the source-drain current and GBL-FET charac- 
teristics are determined by the conductivity of the gated 
section, which, in turn, is determined by the electron den- 
sity and scattering mechanisms including the electron- 
electron scattering mechanism, and by the self-consistent 
electric field directed in the channel plane. In such a 
situation, different hydrodynamical models of the elec- 
tron transport (including the drift-diffusion model) can 
be applied (see, for instance, Refs. 0, H, [H, HI H3 ) . If 
Vth 2 < Vt < Vth.i, considering the potential distribution 



3 



in the direction perpendicular to the GBL plane invoking 
the gradual channel approximation [28l [29 ] and assum- 
ing for simplicity that the thicknesses of the gate layers 
separating the channel and the pertinent gates, Wb and 
Wt, are equal to each other Wb — Wt — W, we obtain 



A = 



(V b + V t ) 



5A(t) 



1 



:6V(t), (2) 



where e = |e| is the electron charge. In the voltage range 
in question, the electron system in the gated section is not 
degenerate. This voltage range as well as the range Vt < 
V th ,2 < Vth.i correspond to the GBL-FET "off-state". 
Similar formulas take place for the barrier height from 
the drain side (with the replacement of A by A + eVd), 
In the cases when V t h,2 < V t h,i < V t or V t < V t h,2 < 

A = -e ± — - — £+, 



1 2-ireW 

SA(t) = --eSV(t)±—--S^{t). (3) 

2 K 

Here + SH^ft) are the electron and hole densities 
in the gated section and k is the dielectric constant of 
the gate layers. In the most interesting case when the 
electron densities in the source and drain sections are 
sufficiently large, so that the electron systems in these 
sections are degenerate. Considering this, the height of 
the barrier A is given by 



A = 



A 



V t (a B /8W) 
'[l + {a B /AW)] 



-eV t 



( aB 
\8W 



V t (a B /8W) + (V t - V b ){d/2W) 
[1 + (a B /4W)] 



(4) 



(5) 



when V t h,2 < V t h.\ < V and V t < V t h,2 < V th ,i, respec- 
tively. Here a B = nh 2 /me 2 is the Bohr radius, d is the 
effective spacing between the graphene layers in the GBL 
which accounts for the screening of the electric field be- 
tween these layers [2(J Hl| . This quantity is smaller than 
the real spacing between the graphene layers in the GBL 
do ~ 0.36 nm. The Bohr radius a B and parametes a can 
be rather different in different materials of the gate lay- 
ers. In the cases of Si0 2 and Hf0 2 (with k ~ 20 [11 1]}) 
gate gate layers, a B ~ 4 nm and a B ~ 20 nm, respec- 
tively. In deriving Eqs. (4) and (5), we have taken into 
account that in real GBL-FETs, (a B /8W) < 1. 
The Fermi energy in the source section is given by 



e F 



[1 + (a B /8W)} 



In 



exp 



a B e 



— V 



8W k B T ) 



1 



(6) 



where k B is the Boltzmann constant and T is the tem- 
perature, so that at sufficiently large back gate voltages, 



e F ~ eV b 



(a B /8W) 
[1 + (a B /8W)} 



, . , a B 
eVb[ 8W 



> k B T. (7) 



Here we have considered that the electron density in the 
source section = n Vb/ineW (the electron density in 
the drain section of the channel is approximately equal 
to = n(V b - V d )/4ireW). Comparing Eqs. (2), (4), 
and (5), one can see that the height of the barrier Ao 
increases with increasing absolute value of the top-gate 
voltage rather slow in the voltage ranges "a" and "c" in 
contrast to its steep increase in the voltage range "b" . 

Since the energy gaps in GBLs E 3iS , E n and En a de- 
pend on the local transverse electric field [20, [2J, [23] , so 
that they are different in different sections of the channel 
depending on the bias voltages: 



E„ 



edV b 



E„ 



ed Q (V b -V t ) 
2W 



E, 



g,d 



edjVb - V d ) 
2W 



(8) 

One can see that at V t < and Vd > 0, one obtains 

Eg > Eg tS > Eg^d- 

The threshold voltages V t h,i and V t h,2 are deter- 
mined by the conditions A = e F and A = s F + 
Eg ,respectively. The latter implies that the Fermi energy 

of holes in the gated section e^ o/e ^ = Ao — e F — E g = 0. 
As a result, the threshold voltages are given by 



V th , 



-Vb[ 1 



w)> Vth < 2 ^~ Vb { 1+ ^ + W 



(9) 

Since one can assume that d <§C W, the threshold voltages 



are close to each other: \Vt 



th,l\ 



< w 



th,2 



with \Vth,2 



V t h,i\ - (2d /W)V b > ^Eg^/e. The values of the energy 
gap in the gated section at the threshold top gate voltages 
are given by 



E,. 



Vt=V th ,i 



Vt=V th - 



W 



2E, 



(10) 



In the following we restrict our consideration by the 
situations when the height of the barrier for electrons 
in the gated section is sufficiently large (so that Ao > 
e F ), which corresponds to the band diagrams shown in 
Figs. 2(b) and 2(c). 



III. BOLTZMANN KINETIC EQUATION AND 
ITS SOLUTIONS 

The quasi-classical Boltzmann kinetic equation govern- 
ing the electron distribution function / p = / p (x, t) in the 
section of the channel covered by the top gate (gated sec- 
tion) can be presented as 



df P , df p 



dt 



c). 



J d 2 qw(q)(f p+cl - f p )8(e p+cl 



(11) 

Here, taking into account that the electron (and hole) 
dispersion relation at the energies close to the bottom of 
the conduction band is virtually parabolic with the effec- 
tive mass m (m ~ 4 x 10~ 29 g), for the energy of electron 
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with momentum p = (p Xl p y ) we put e p = p 2 /2m = e, 
v x = Px/ m = pcosO/m, where cos# — p x /p (the x-axis 
and the y-axis are directed in the GBL plane) and w(q) 
is the probability of the electron scattering on disorder 
and acoustic phonons with the variation of the electron 
momentum by quantity q. 

The density of the electron (thermionic) current, J = 
J(x,t), in the gated section of the channel (per unit 
length in the y-direction) can be calculated using the 
following formula: 



J = 



4e 



(2-Kh) 2 



d pv x f. 



(12) 



where h is the reduced Planck constant. 

Disregarding the electron-electron collisions in the 
gated section of the channel (due a low electron den- 
sity in this section in contrast to the source and drain 
sections where the electron-electron collisions are essen- 
tial), we consider two limiting cases: ballistic transport of 
electrons across the gated section and strongly collisional 
electron transport. 



A. Ballistic electron transport 



Q{p x ) in Eqs. (16) and (17) reflects the fact that there are 
no electrons propagating backwards due to the absence 
of the electron scattering in the gated section. 

Using Eqs. (12), (14), and (15), we arrive at the fol- 
lowing formulas for the dc and ac components, Jo and 
S Ju , of the current at the drain edge of the gated section 
(i.e., at x = L t ) 



Jo 



A 



'0 > 



SJu, 

Jn 



6A U 
k B T 



B-L J JO 



Here 



„ i 1/4 ( 

2^= exp i— = 



(16) 



Jo 



C 



2v^ 



iS 



2V? 



If SV(t) = SV^e-™*, where 5V U < \V t \ and u are 
the amplitude and frequency of the ac signal, then the 
electron distribution function can be searched as f p = 



F a + SF u {x) e- luJt and A = A + <5A W e 



Assuming 



that eVd 3> ksT and solving Eq. (11) with the boundary 
conditions 



exp 



e F - A(t) - e 



knT 



p x >0,x=0 

where L t is the top gate length, we obtain 
' e F - A - e 



0. 



p x <0,x=L t 



F ~ exp 



k B T 



(13) 
(14) 



exp 



e F - A - e 
k R T 



2e cos 9 



5A^ 
k B T 



e(p B ). 



(15) 



Here, Q(p x ) is the unity step function. The first bound- 
ary condition given by Eq. (14) corresponds to quasi- 
equilibrium electron distribution in the source section of 
the channel and the injection of electrons with the kinetic 
energy exceeding the barrier height A(i) from the source 
section to the gated section (at x = 0). The injection 
of electrons from the drain source to the gated section 
(at x — Lt) is negleted due to eVd ^> ksT; this inequal- 
ity leads to rather high barrier near the drain edge of 
the gated section. The presence of the unity step fuction 



where t = Lty / m/2kBT is the effective ballistic transit 
time across the gated section of electrons with the ther- 
mal velocity v F — ^j2ksT /m and C(x) and S(x) are 
Frenel's cosine and sine functions. At lot 3> 1, 



At lot <C 1, tends to unity. 



B. Collisional electron transport 

In the case of strongly collisional electron transport, 
the distribution function in the gated section is close to 
isotropic and it can be searched in the form 



f p = F + g cos 0. 



(18) 



Here F — F(e, x, t) is the symmetrical part of the elec- 
tron distribution function (which is generally not the 
equilibrium function). The second term in Eq. (18) 
presents the asymmetric part of the distribution func- 
tion with g = g(e,x,t). Similar approach was used for 
the calculation of characteristics of heterojunction bipo- 
lar transistors [32j (see also Refs. PHHH). As a result, 
after the averaging of Eq. (11) over the angle 9, one can 
arrive at the following coupled equations: 



dF 
~dt 



e dg 
2m dx 1 



dgF 
dt 



vg = 



m ox 



■5 



Here in the case of w(q) — w — const, which corre- 
sponds to the scattering of electrons on short-range de- 
fects, v — mw/2. Equations (19) are reduced to the 
following equation for function F: 



d_fdF e 



dt \ dt 



vF 



e d 2 F 
mdx 2 ' 



(20) 



In the most interesting case when eVd 3> ^T, the 
boundary conditions for Eq. (20) at x = and x = L t , 
can be adopted in the following form: 



= cxp 



x=0 



e F - A(t) - e 



k B T 



0. (21) 



x=L t 



cosh[a ul (x — L t )\ ( 5 A 



smh.(a u L t ) 



(27) 



After that, using Eqs. (12), (26), and (27), we arrive at 
the following formulas for J° and SJ^ (at x = L t ): 



Here 



2{k B T) 2 fe F -A 
Jo = e — —r- cxp 



Jo 



irh 2 L t v 



5A a 



k B T 



— J i 



k B T 



pOO 

Jo 



(28) 



(29) 



The boundary condition under consideration imply that 
at x = there is the electron injection from the source 
section of the channel, whereas at x = L t an effective 
extraction of the electrons into the drain section occurs 
due to a strong pulling dc electric field. Due to a strong 
electron scattering a significant portion of the injected 
electrons returns back to the source section. 

Setting as above SV(t) = SV^ e~ iut and, hence, F = 



5F U 



and g = F Q + Sg^ e~ 



we obtain 



d 2 F 
dx 2 



= 0, 



d 2 SF LU muj{uLi + iv p ) 



dx 2 



and arrive at 
SF^ = cxp 



exp 



S F 



k B T 



5F„ = 0, 



x 

1 



(22) 



(23) 



(24) 



e F - A 



k B T 



sinh[a w (x - L t )\ (5 A 



sinh(a w L t ) 



knT 



(25) 

where = y/mLd(uj + iu)/e. Considering Eqs. (19), 
(24), and (25) , we obtain 



l TV 



sinh ^[2u(uj + iv)r 2 /£] Vw + if' 



2lj 



According to Eq. (28), Jq oc 1/L t v. One needs to stress 
that the collisional case under consideration corresponds 
actually to > 1. In the frequency range ui <C v/t 2 , 
~ V?- At u> > v/t 2 , one obtains 



^exp 



(1 + i)\Jllit 2 v 

71 . 



IV. GBL-FET DC TRANSCONDUCTANCE 

Equations (16) and (28) provide the dependences of 
the source-drain dc current J as a function of the de- 
vice structural parameters, temperature, and back- and 
top-gate voltages for GBL-FETs with ballistic and colli- 
sional electron transport, respectively (in the limit eVd ^> 
k B T). Using Eq. (16), one can find the dc transconduc- 
tance G = (dJ /dV t )\ Vb of a GBL-FET with the ballistic 
electron transport: 



B _ 2 \/2mk B T 



exp 



£F - A 
k B T 



2k B T 



(30) 



when Vth.2 < V t < V th ,i, and 



90 



1 2edFo 
v V m dx 



1 2s 

— exp 
i>L t \m 



e F - A 



= — zcxp 



k B T 



2ed5F LU 



(lo + iv)\ m dx 



£ F - A - 
k B T 



2uj 



(lu + iv) 



(26) 



,y/2mk B T 



£F - A 



when V t < V t h.2 < Vth.i- Here 



eJl_ 
2k B T 



(31) 



(32) 



Similarly, using Eq. (28), we obtain the following for- 
mulas for the GBL-FET transconductance in the case of 
collisional electron transport: 



; k B T 



exp 



k B T 



eJ§_ 
2k B T 



(33) 
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when Vth.2 < V t < V th ,i and 



irh 2 L*v 



exp 



n 



2k B T 



K (34) 



when V t < V t h,2 < Vth,\- 

As follows from the comparison of Eq. (30) with 
Eq. (31) and Eq. (33) with Eq. (34), the GBL-FET dc 
transconductance in the top gate voltage range V t h,2 < 
V t < Vthi (in the range "b") might be much larger than 
that when V t < V t h,2 < V t h,i (in the range "c") since 
as <C 8W. This is due to relatively slow increase in Ao 
with increasing \Vt\ when the hole density in the gated 
section becomes essential. 

The voltage dependences of the dc transconductance 
can be obtained using Eqs. (30), (31), (33), and (34) and 
invoking Eqs. (2), (4), and (6). In particular, in a rather 
narrow voltage range Vth,2 <Vt< Vth,i, one obtains 



G C = G B K exp 

VT 



e(V t - V thtl ) 



2k B T 



(35) 




COX 



FIG. 3: The ac transconductance (normalized by Gq ) versus 
tor calculated for GBL-FETs with ballistic (v = 0) and col- 
lisional electron transport (vt = 2,4, and 6) at the top gate 
voltage in range "b", i.e., Vth,i <Vt < Vth,l- 



At sufficiently large absolute values of the top-gate volt- 
age when Vt < Vth,2 < Vth,i, the transconductance vs 
voltage dependence is given by 



n C n B ( a B \ 

Go =— G„ oc^— jexp 



e(V t ~ V b ) ( d 
k B T \2W 



x exp 



e(Vt - Vthg) ( ob_ 
\8W 



k B T 



(36) 



V. GBL-FET AC TRANSCONDUCTANCE 

According to the Shockley-Ramo theorem [34, 35], the 
source-drain ac current is equal to the ac current in- 
duced in the highly conducting quasi-neutral portion of 
the drain section of the channel and in the drain contact 
by the electrons injected from the gated section. This 
current is determined by the injected ac current given 
by Eq. (16) or Eq. (28) as well as the electron transit- 
time effects in the depleted portion of the drain section. 
However, if tOTd -C 1, where the Ta is the electron transit 
time in depleted region in question, the induced ac cur- 
rent is very close to the injected ac current [H, HE H3- 
Since, at moderate drain voltages, the length of depleted 
portion of the drain section Ld can usually be shorter 
than the top gate length L t , in the most practical range 
of the signal frequencies w < r" 1 , one can assume that 
Td <C r and, hence, uTd <C 1. Considering this and us- 
ing Eqs. (17) and (28), the GBL-FET ac transconduc- 
tance G w = (dS J w /dSV w )\v b at different electron trans- 
port conditions can be presented as 



G. 



85A U 



k R T\ dSVu 



/ d^e-^O, (37) 
Jo 



G, 



c _ J o 



dSA u 



k R T\ 85 Vu, 



/ d^e^HUO, (38) 
Jo 



respectively. 

In the range of gate voltages V t h,2 < V t < V t h,i (range 
"b"), Eq. (2) yields 



d5K 



v d 



ddVu 

In this case, Eqs. (37) and (38) result in 
eJa 8 



2k B T 



G^ = 



eJ 

2k B T 



Jo 



(39) 



(40) 



(41) 



As follows from Eqs. (40) and (41), the characteristic 
frequencies of the ac transconductance roll-off are 1/r 
and v/t 2 in the case of the ballistic and collisional elec- 
tron transport, respectively, i.e., the inverse times of the 
ballistic and diffusive transit across the gated section of 
the channel. Indeed, the quantity v/t 2 can be presented 
as D/Lf, where D is the electron diffusion coefficient. 

The situation becomes more complex in the range of 
the top gate bias voltages V t < V t h,2 < V t h,i (range "c"). 
As follows from Eq. (4) in this voltage range, the quan- 
tity 5 is determined not only by the ac voltage 6 VL> 
but also by the ac component of the hole density in the 
gated section <$£+. Moreover, at sufficiently high signal 
frequencies, the hole system in the gated section can not 
manage to follow the variation of the ac voltage. Taking 
into account the dynamic response of the hole system 
(see Appendix A), instead of Eq. (39) one can obtain 
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FIG. 4: The ac transconductance (normalized by Gjf ) versus 
lot calculated for GBL-FETs with collisional electron trans- 
port (vt — 2) at the top voltage in range "c" ( Vt < V t h,2 < 
V t h,i) at T r /r = 50. 



dSA^ 



dSV u 

Here (see Appendix A) 



(42) 



v d 



{(«) + _L_r' 

\\4Wj 1-iWTr) 



1 



lLOT r 



1 — ILJTrlZo 



(43) 



where r r is the time of the gated section recharging asso- 
ciated with changing of the hole density due to the tun- 
neling or/and generation- recombination processes. Gen- 
erally, r r depends on the top gate length L t . 

Accounting for Eq. (42), we arrive at the following 
formulas for the GBL-FET ac transconductance when 
V t < V th ,2 < V th x- 



G c , = 



eJ Q B 
2k B T 



poo 

Jo 



e-^ w (0, (44) 



^o „ 
2k B T u 



Jo 



(45) 



If uj 3> TjT 1 , one obtains 1Z U ~ 1 and the ac transconduc- 
tances in both "b" and "c" ranges of the top gate voltage 
are close to each other (compare Eqs. (40) and (41) with 
Eqs. (42) and (44)). However, at low signal frequencies 
( uj 3> TjT 1 ), the ac transconductance given by Eq. (44) 
or (45) for the voltage range "c" are markedly smaller 
than those given by Eqs. (40) or (41) valid in the voltage 
range "b". 



VI. ANALYSIS OF THE RESULTS AND 
DISCUSSION 

Comparing Gq and Gq given by Eqs. (30) and (35), 
we obtain Gq / Gq = ^/tt/ vt. This implies that the above 
ratio markedly decreases with increasing collision fre- 
quency (with decreasing electron mobility) and the top 
gate length, i.e., with the departure from the ballistic 
transport. 

As shown above, the dc current steeply drops in a nar- 
row top-gate voltage range V t h.2 <V t < V t h,\- Indcd, the 
ratio J B \v t =v th , 2 /Jo\v t =v thA ^ exp[-(edoV b /Wk B T)]. 
Setting W = 5 nm, T = 300 K, and V b = 1 - 2 V, we find 
Jo\v t =v th ,JJ B \v t =v thA - 3 x 10- 3 -6 x 10- 2 . The esti- 
mate of the dc current at Vt = Vth.2 — —Vb (which might 
be interesting for the GBL-FET applications in digital 
large scale circuits) with W = 5 nm, T = 300 K, and 
V b = 1 - 2 V yields J S ~ 1 x 10~ 3 - 2 x 10~ 2 A/cm. At 
T = 77 K and V b = 1 V, one obtains J S ~ 7x 10" 7 A/cm. 
In the case of a GBL-FET with the width H = 1 fim, the 
latter corresponds to the characteristic value of the "off- 
current" Jq H ~ 70 pA. Similar values can be obtained 
at T = 300 K when -V t = V b ~ 4.6 V. 

As follows from Eq. (35), the GBL dc transconductance 
in the range of the top gate voltages V t h,2 < Vt < V t h,\ 
is particularly large when V t < V t h,i- This is due to 
a sharp voltage- sensitivity of the dc current and its 
relatively high values at such voltages. Indeed, using 
Eq. (28) at T = 300 K and V t < V th ,i, we obtain 
G B < 2500 mS/mm. In GBL-FETs with W t < W b , 
the dc transconductance can be even larger. 

The pre-exponential factor in the right-hand side of 
Eq. (36) is proportional to a small parameter (clb/8W). 
The argument of the exponential function in this equa- 
tion comprises small parameters (cib/&W) and (do/2W). 
This implies that the dc transconductance in the volt- 
age range V t < V t h,2 < V t h,i described by Eq. (36) is 
relatively small and is a faily weak function of the top- 
gate voltage. As follows from Eqs. (35) and (36), the 
ratio of the dc transconductance at Vt < Vth.2 to that 
a t V t < V t h,i is equal approximately to the following 
small value:' {a B /4W)exp[edo(V t h,2 - V b )/2Wk B T} = 
(aB/4:W)exp(—2eEg^do/dkBT); it is smaller then the 
ratio of the dc currents by parameter (a_e/4H r ). 

Figure 3 shows the ac transconductance \G U \ normal- 
ized by the dc transconductance at the ballistic trans- 
port Gq as a function of the normalized signal frequency 
lot calculated for GBL-FETs with both ballistic and col- 
lisional electron transport. It is assumed that the top 
gate voltage is in the range V t h,2 < Vt < V t h,i- The 
inset in Fig. 3 shows the dependence of the normalized 
threshold frequency LOtT on vt oc v L t . The threshold 
frequency is defined as that at which \G u \/Go = l/y/2. 
One can see that \G U \ pronouncedly decreases with in- 
creasing collision frequency. However, as seen from the 
inset in Fig. 3, the decrease in ui t with increasing v is 
markedly slower: the ratio of Lo t at vt — 2 and u) t at 
vt = 6 (i.e., three times larger) is approximately equal 
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to 1.42. Setting L t = 100 - 500 nm and T = 300 K, 
for the threshold frequency f t = u t /2n at the ballis- 
tic transport we obtain f B ~ 0.485 - 0.97 THz. To 
realize the near ballistic regime of the electron trans- 
port (yr -C 1) in GBL-FETs with such gate lengths, 
the electron mobility /i > (1 — 5)xl0 4 cm 2 /Vs is re- 
quired. The possibility of the latter mobilities at room 
temperatures was discussed recently (sec, for instance, 
Ref. 0). At a shorter top gate, L t = 75 nm, one obtains 
f B ~ 1.29 THz. In the case of a GBL-FET with rela- 
tively long top gate and moderate mobility (L t = 500 nm 
and n = 2 x 10 4 cm 2 /Vs) when the effect of scattering 
is strong (y ~ 2 x 10 12 s _1 and vt ~ 2.7), we obtain 
ff ~ 94 GHz. 

In Fig. 4, the similar dependences calculated for a 
GBL-FET with collisional electron transport at Vth.2 < 
Vt < Vth.i and at V t < V t h,2 < V t h,i are demonstrated. 
Since in the latter top-gate voltage range the electron 
scattering on holes accumulated in the gated section can 
be strong, so that the realization of the ballistic trans- 
port at such top-gate voltages might be problematic, only 
the dependences corresponding to the collisional electron 
transport are shown. As follows from Eqs. (44) and (45) 
and seen from Fig. 4, the ac transconductance at the top 
gate voltages V t < V t h,2 < Vth,i is fairly small at low 
frequencies u> < r" 1 being close to the dc transconduc- 
tance (due to a smallness of parameter 1Z ~ as/4W), 
whereas it becomes much larger in the intermediate fre- 
quency range t^ 1 <C lo < t _1 (or t^ 1 <C lo < vt~ 2 ). 
This is due to the effect of holes in the gated section. 
Owing to this effect, the low frequency noises can be ef- 
fectively suppressed. 

Using the above results for GBL-FETs and those ob- 
tained previously [l9[ for GNR-FETs, one can compare 
the GBL and GNR characteristics. In particular, consid- 
ering the expressions for J B found for GBL-FETs and 
GNR-FETs, we obtain J^/J R ' GNR = 1/2. For the case 
of collisional transport, one obtains J G / J G ' GNR ~ 1. As 
a result, the GBL-FET and GNR-FET dc transconduc- 
tances are close to each other. The ratio of the GBL-FET 
and GNR-FET ac transconductances at high signal fre- 
quencies is G B /G B ' GNR oc l/V^, i.e., the GBL-FET 
ac transconductance falls more steeply with increasing lo 
than the GNR-FET ac transconductance. 

The GBL-FET dc and ac characteristics obtained are 
valid if the interband tunneling source-drain current 
through the n + — p and p — n + junctions beneath the 
edges of the top gate are small in comparison with the 
thermionic current created by the electrons overcoming 
the potential barrier in the gated section. Such a tun- 
neling current can be essential in the voltage range "c" 
(V t < V t h,2 < V t h,x) depending on the energy gap near 
the top gates and the length of the n + — p and p — n + 
junctions in question. This implies that there is a lim- 
itation when the top-gate voltage is not too high Vt in 
comparison with the threshold voltage Vth.2 (he., when 
Vt < Vthfi), so that the calculated characteristics corre- 
spond to the most interesting voltage range where the ac 



transconductance can be rather large. 

VII. CONCLUSIONS 

We have presented an analytical device model for a 
GBL-FET. Using this model, we have calculated the 
GBL-FET dc and ac characteristics and shown that: 

(1) The dependence of the dc current on the top gate 
voltage is characterized by the existence of three voltage 
ranges, corresponding to (a) the population of the gated 
section by electrons, (b) the depletion of this section, and 
(c) its essential filling with holes, and determined by the 
top-gate threshold voltages V t h,i and V t h,2- 

(2) The ac current is most sensitive to the top-gate volt- 
age Vt and the dc and ac transconductances are large 
when Vth.2 < Vt < Vth.2- 

(3) The electron scattering in the gated section results 
in a marked reduction in the dc and ac transconduc- 
tances. However, the threshold frequency corresponding 
to \G u \/Go = l/v 7 ^ decreases with increasing collision 
frequency relatively smoothly. 

(4) The transient recharging of the gated section by holes 
(at Vf < Vth.2) leads to a nonmonotonic frequency de- 
pendence of the ac transconductance with a pronounced 
maxima in the range of fairly high frequencies. This ef- 
fect might be used for the optimization of GBL-FETs 
with reduced sensitivity to low frequency noises. 

(5) The fabrication of GBLs with high electron mobility 
at elevated temperatures opens up the prospects of re- 
alization of terahertz GBL-FETs with ballistic electron 
transport operating at room temperatures and surpass- 
ing FETs based on A 3 B 5 compounds. 

Acknowledgments 

The authors are grateful to Professor S. Brazovskii 
(University Paris-Sud) and Professor E. Sano (Hokkaido 
University) for fruitful discussions and valuable informa- 
tion. The work was supported by the Japan Science and 
Technology Agency, CREST, Japan. One of the authors 
(N.K.) acknowledges the support by the INTAS grant 
7972, and by the ANR program in France (the project 
BLAN07-3-192276). 

Appendix A. Dynamic response of the hole system 
in the gated section 

At sufficiently large values |V t |, the gated section is 
essentially populated with the holes. As follows from 
eq. (4), 

SA u = --e5V u «$£+. (Al) 

i K 

The ac component of the hole density <5£+ obeys the con- 
tinuity equation which can be presented in the following 
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form: 



:SS+ = SG U 



8J^\ X=Q -5JT,\,= 



u> \x=L t 



eL t 



(A2) 



Here SG^ is the variation of the generation pf holes 
in the gated section (associated, say, with the generation 
of holes by the thermal radiation [381 ]) and 5 J^\ x= o and 
8 J^\ x =L t are the interband tunneling ac currents near 
the source and drain edges of the top gate, respectively. 
For normal operation of GBL-FETs, these tunneling cur- 
rent should relatively small. This is achieved in GNL- 
FETs by proper choice of the energy gap in the different 
sections of the channel [E g s , E g „ and E g( i which should 
be not too small. The terms in thr right-hand side of 
eq. (A2) can be presented as 



5G U = K g (SA u - Ss+), 



5.1 



ui \x=0 



61 



ui \x=L t 



2{5A UJ - 6e+) 



eR t 



(A3) 



(A4) 



Here SeJ, ~ (7rfi, 2 /2m)<51]+, K g = 2m/irfi 2 Tg, where r g is 
the characteristic time of the generation-recombination 
processes in the gated section, R t is the tunneling resis- 
tance of the p-n junctions induced by the negative top 
gate voltage near the edges of the top gate: lnR t oc 
El ,2 /E\\ oc {V b ~Vtf/ 2 /V t . Fromeqs. (Al) - (A4), taking 
into account the limit u — > (see eqs. (5) and (30)), we 
find 



S A u 



■\sv u 



( a B 
x 



Ob_\ 

AW J 



1 



(a B \ 1 

W T RC 



1 (a B \ 1 iu 

Tg \AW J T RC 



i 



(A5) 



where trc — RtCt is the time of the gated section 
recharging by the tunneling currents, R t is the tunnel- 
ing resistance of the p-n junctions induced by the nega- 
tive top gate voltage near the edges of the top gate, the 
quantity r r = trct 9 /[trc + (as / 4W)r g ] is the charac- 
teristic time of the gated section recharging by holes, and 
Ct = KL t /2irW is the capacitance of the gated section. 
At oj — > 0, eq. (A5) leads to 



where 



1 



-eSV^TZo, 



dSVuj 



-TZ , (A6) 



TZ 



{a B /W) 



(a B /4W), 



1 + (a B /4W) 

i.e., coincides with the value given by eq. (28). When 
u) ^> TjT 1 , eq. (A5) yields 



5A U 



:eSV a , 



dSA u 



dSVuj 



(A7) 
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